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Abstract: The hadronic decays rf — > tjtttt are studied in the frameworks of large- Nc 
Chiral Perturbation Theory, at lowest and next-to- leading orders, and Resonance Chiral 
Theory in the leading 1/Nc approximation. Higher order effects such as irir final state in- 
teractions are taken into account through a detailed unitarization procedure. The inclusion 
of finite-width effects in the case of RChT is also discussed. The Dalitz plot distribution 
and the differential branching ratio are computed in both approaches. The predicted Dalitz 
plot parameters obtained from the different treatments are compared with the most recent 
measured values. We find that the rf — > rjinr branching ratios are easily understood, while 
the Dalitz plot parameters require the inclusion of tttt loops in order to achieve a reasonable 
agreement. Our final predictions agree with the experimental measurements. We hope our 
results to be of relevance for present and future experimental analyses of these decays. 
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1. Introduction 

Chiral Perturbation Theory (ChPT) Q is the low-energy effective theory of Quantum 
Chromodynamics (QCD). It is built from the global SU(3)l x SU(3)r symmetry exhibited 
by QCD in the chiral limit, and described in terms of an octet of pseudoscalar bosons 
appearing in the theory, as a result of the spontaneous breaking of this symmetry down 
to SU(3)y. The eight Goldstone bosons are identified with the lightest hadronic states 
(n, K, rj) and their small masses are generated from the quark mass term, which explicitly 
breaks the global symmetry of QCD. The ChPT Lagrangian is organized in terms of an 
increasing number of powers of momentum and quark masses. Below the resonance region 
{E < M p ), the interactions of the (jr,K,r]) particles are systematically analyzed and easily 
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understood within this framework. The enormous success in the description of these low- 
energy interactions makes of ChPT a powerful theoretical tool Q. In this formalism, 
however, the pseudoscalar singlet rji is not explicitly included since the U(1)a anomaly 
prevents this state from becoming a ninth Goldstone boson. Therefore, processes involving 
the rji as an external degree of freedom are not accounted for in ChPT 1 . This drawback is 
remedied in the large- Nq limit, where the effects of the axial anomaly are absent and the 
symmetry is enlarged to U(3)l x U(3)r. Then, a simultaneous expansion in p 2 , m q and 
X/Nq is possible and the interactions among the (ir, K, rj, rf) mesons can be described with 
a Lagrangian. Whether this extended framework, named large- Nc ChPT after the work 
of Kaiser and Leutwyler Q, is well established is still under discussion due to the large 
physical mass of the if. In addition, large- Nc ChPT does not include resonances as external 
states. Their masses are supposed to be bigger than the center-of-mass total energy of a 
given process involving the pseudoscalar mesons and consequently they are integrated out. 
The effects of these resonances are then virtual and encoded in the low-energy constants 
of the chiral Lagrangian. However, when the energy of the process is of the order of the 
resonance mass, the perturbative expansion of ChPT stops being valid and the resonant 
effects must be taken into account explicitly. This is considered in Resonance Chiral Theory 
(RChT) H, where the interactions of the pseudoscalar mesons are supplemented with new 
interactions among these and nonets of vectors, axials, scalars, representing the p, a±, a, 
etc., in a minimal way. 

The decays rj' —> rjirn are interesting for several reasons. First, due to the quantum 
numbers of the pseudoscalar mesons involved, any resonance involved in the decay must 
be mainly of scalar nature. G-parity prevents vectors from contributing. Therefore, this 
decay is specially suitable for the analysis of the properties of the /o(600) (or a) resonance, 
even though the ao(980) is also present and, in fact, dominant. Second, the presence of 
77 and rf in this reaction is ideal for studying the mixing properties of these two mesons. 
Third, and more general, this decay allows to test ChPT and its possible extensions such as 
large- Nc ChPT and RChT. In view of all that, precision measurements on rj and rf would 
be very helpful and would provide useful information on our understanding of low-energy 
QCD. In particular, their decays are very useful for studying symmetries and symmetry 
breakings in QCD. The simultaneous treatment of both the 77 and rf imposes constraints 
on theoretical approaches which are tighter than those considering the rj alone. 

Accordingly, there is at present an intense activity studying these processes. Recently, 
the GAMS-47T and VES Collaborations have measured the related Dalitz plot parameters 
(GAMS-47T for the rf — > 7?7r°7r° channel || an d VES for the rf — > r/7r + 7r~ one Q) comple- 
menting older results reported by an early GAMS Coll. |7j and CLEO Q. In the isospin 
limit the values of the Dalitz plot parameters should be the same but the experimental 

1 Strictly speaking, ChPT takes into account the virtual effects of the pseudoscalar singlet by means of 
the low-energy constants. 



- 2 - 



measurements show some discrepancies among them. Therefore, there is a clear need for 
improvement in experimental precision and new measurements of these parameters are fore- 
seen at KLOE-2, Crystal Ball, Crystal Barrel and maybe WASA with improved statistics 
and a good understanding of the systematical errors. At KLOE the rf is produced via the 
process e + e~ — > <!> followed by <I> — > 7/7. For the rf —> rjir + ir~ decay channel around 21K 
events are now on tape. The background is seen to be very low and these data can easily 
be used to determine the Dalitz plot parameters. The increased luminosity at KLOE-2 (a 
factor 3) will open new possibilities for rf studies. For instance, a Monte Carlo simulation of 
the rf — > r]Tr + ir~ process shows that the detector has a good sensitivity to the a meson ||. 
At Crystal Ball, the rj' is photo-produced, jp — > rj'p, and can be identified, e.g., from its 
decay rf — > r/7r°7r° — > 67. About 10K events were expected during 2009 for the rf — > ry7r°7r° 
process. At WASA, large samples of rf will be produced in the reaction pp — > pprf as soon 
as the new detector is developed. Finally, about 27 million and 12 million decay events 
could be detected at BES-III each year for rf — > r]ir + ir~ and rj' — > rjir 7r°, respectively 
1 1 C ] . On the theory side, the rj' — > rjirir decays have been studied within an effective chiral 



Lagrangian approach in which the lowest lying scalar mesons are combined into a possible 
nonet more recently, within the framework of U(3) chiral effective field theory in 

combination with a relativistic coupled-channels approach |12[ |. Older analyses based on 
chiral symmetric frameworks include contact terms |l3|, [TIL [l5|, [lTj or contact terms 



plus scalar meson exchanges [18, IS]. 



In this work we study the decays rj' — > rjixis in two well established chiral frameworks, 
large- Nq ChPT and RChT. Since the charged decay is the most prominent channel of 
the rj' we perform the whole analysis for this decay. Working in the isospin limit, the 
predictions for the neutral channel will be the same at the level of amplitudes as for the 
charged channel, and one half at the level of the branching ratios. Preliminary results were 
presented in Refs. p0|, 



The paper is divided as follows. In Sec. 0, we introduce the Dalitz plot parameterization 
and discuss the different values of the related parameters. We also report on the latest 
experimental values of the charged and neutral branching ratios together with their ratio. 
Sec. H is devoted to the analysis of the rf — > rjir^-K^ process in the framework of large- 
Nq ChPT, first computing its predictions at leading and next-to-leading order and then 
considering the corrections due to the unitarization of the tttt channel. In Sec. |], the 
decay rj' — > r/7r + 7r~ is analyzed in the framework of RChT starting with the leading order 
prediction and later calculating the subleading corrections due to resonance widths and 
the unitarization of the tttt channel. Finally, we discuss the results obtained and present 
our conclusions in Sec. ||. The techniques of partial-wave projection and unitarization of 
final state interactions (FSI) used in this work along with a description of rj-rf mixing for 
states and decay constants are included in the appendices for completeness. 



-3- 



2. Experimental status of Dalitz plot parameters and branching ratios 



The Dalitz plot distribution for the charged decay channel is described by the following 
two variables: ^ 

X --Q(T n+ -T 7r -), Y — — - — - 1 , (2.1) 

where T n ± „ denote the kinetic energies of mesons in the rj 1 rest frame, 

(m v > - m v ) 2 -s (m v > - m w ) 2 - t {m n > - m^) 2 -u 

+ v = o 5 i tt+ = 7, > 1 n- = 7, j \ 2 - 2 ) 

Zm^i Zmw Irrirf 

and Q = T V + T n + + T n - = — — 2m n . The Mandelstam variables s = (p 7T + +p 7T -) 2 = 
m 2 ^, t = (prj/ — p n +) 2 = m 2 and u = (p v > — p n -) 2 have been employed here, which obey 
the relation s + t + u = m 2 , + m 2 + 2m 2 . 

The squared absolute values of both decay amplitudes (charged and neutral) are ex- 
panded around the center of the corresponding Dalitz plot in order to obtain the Dalitz 
slope parameters 2 ||: 

|^(X,y)| 2 = \M\ 2 [l + aY + bY 2 + cX + dX 2 } , (2.3) 

where a, b, c and d are real parameters and \M\ 2 is a normalization factor. For the charged 
channel odd terms in X are forbidden due to charge conjugation symmetry, while for the 
neutral c = from symmetry of the wave function. The Dalitz plot parameters may be 
different for charged and neutral decay channels. However, in the isospin limit they should 
be the same. A second parametrization is the linear one 

\A(X,Y)\ 2 oc \l + aY\ 2 + cX + dX 2 , (2.4) 

where a is a complex parameter. Comparison with the general fit gives a = 2Re(a) and 
b = Re 2 (a) + Im 2 (a). The two parametrization are equivalent if b > a 2 /4. 



Parameter 


Exp. [rj' -> riTT°TT°] 
GAMS-4tt [p| 


Th. [rf -5- 7?7r°7r°] 
Borasoy & Nissler jl^] 


Exp. [■>]' — ¥ r/ir+TT ] 
VES 


Th. [rf -> rin+TT-] 
Borasoy & Nissler |l2[ 


a 


-0.066 ±0.016 ±0.003 


-0.127 ± 0.009 


-0.127 ±0.016 ±0.008 


-0.116 ±0.011 


b 


-0.063 ±0.028 ±0.004 


-0.049 ± 0.036 


-0.106 ±0.028 ±0.014 


-0.042 ± 0.034 


c 


-0.107 ±0.096 ±0.003 




±0.015 ±0.011 ±0.014 




d 


+0.018 ±0.078 ±0.006 


±0.011 ± 0.021 


-0.082 ±0.017 ±0.008 


+0.010 ±0.019 



Table 1: Dalitz slope parameters (experiment and theory) for rf — > r]ir ir (second and third 
columns) and rf — > r)it + Tr~ (fourth and fifth columns), respectively. 

The latest available experimental information on the Dalitz slope parameters is sum- 
marized in Table |l]. The analysis by the GAMS-4-7T Collaboration is based on approxi- 
mately 15000 events ||. Note that b is negative here and the fit is thus not compatible 

2 The parameterization in Eq. ( ^.3| ) had been proposed in Ref. (22) with an extra term eXY . The analysis 
of Ref. M included this term in their fits and found that parameter e is consistent with zero. 
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with the linear fit of Eq. fl2.4|) . If a fit is done with the linear parameterization one gets 
a = —0.042 ± 0.008, which is in agreement with an early measurement from GAMS based 
on 5400 events that gave a = —0.058 ± 0.013 (assuming Im(a) = and c = 0) J?j]. Both 
analyses are for the rf — > ^7r°7r° channel. Regarding the rj' — > f]TT + ir^ decay, the re- 
sult from CLEO based on 6700 events yields a = —0.021 ± 0.025 (assuming Im(a) = 0, 
c = and d = 0) using the same linear fit as for the neutral decay channel [||. The 
VES analysis is based on roughly 13600 events obtained from the charge-exchange reaction 
ir~p — > rj'n and 6500 events from the diffractive-like production reaction tt~N — > r]'ir~N 
Q. A fit using the combined data sets from VES gave the parameter values shown in 
Table g. A gain, b is found to be negative and thus incompatible with a linear fit, while c 
is consistent with zero. The previous work supersedes a first study |]24f , where the r]ir + Tr~ 
Dalitz plot has been investigated with a sample of approximately 7000 events obtaining 
a = -0.072 ± 0.012 ± 0.006 —this is the real part of a while Im(a) = 0.0 ± 0.1 ± 0.0— in 
the linear and a = —0.120 ±0.027 ±0.015 in the general parameterization. The C-violation 
parameter was compatible with zero, c = 0.021 ± 0.024. In average, a = —0.059 ± 0.011 



f23| |. Finally, a measurement of a not included in the average is reported in Ref. [25], where 
a = —0.08 ± 0.03 (assuming Im(a) = and c = 0) with about 1400 events. 

For the general parameterization, we see from Table | that there is some tension in all 



the parameters. Notice, however, that the theory model results of Ref. 12] for the charged 



and neutral decay channels are compatible among themselves. While there is agreement in 



a between the VES fitted value and the result of Ref. [12| there is not such when compared 
to the GAMS-47T value. The same happens to d but this time the GAMS-4-7T and the theory 
model results agree, although they are in conflict with the VES reported value. The case 
of b is less severe. The VES value is compatible with the others only at 2a. Finally, there 
is a variation between the measured values of c. Nevertheless, the large statistical error of 
the GAMS-47T result makes this statement not conclusive. 

In this work we want to point out the possibility and maybe the need of extending the 
Dalitz plot parameterization up to higher orders in a systematic way. Hence, we perform 
a general power expansion in Y and X 2 (odd powers of X are forbidden) and consider the 
following extended parameterization 

\A(X, Y)\ 2 = \M\ 2 [l + (aY + dX 2 ) + (bY 2 + k 21 X 2 Y + n m X^) + ■■■ . (2.5) 

Higher terms of the form K2 m ,nX 2m Y n with m + n > 3 are not considered, since they 
are beyond the precision of current and forthcoming experiments. The parameters a, b 
and d have been left with the original nomenclature in order to better compare to former 
analyses. 

Finally, we present the latest experimental values for the charged and neutral rj' — > ijirn 
branching ratios. The CLEO collaboration has recently measured both channels from the 
analysis of J ftp -)■ 7/7 events |2(| and obtained B v <_> V7r+7r ~ = 0.424 ± 0.011 ± 0.004 and 
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0„/_ M p r o 7r o = 0.235 ± 0.013 ± 0.004, respectively. However, we prefer to use the fitted values 
appearing in the most recent version of the Review of Particle Properties [p3fl , that is, 
B n i^ V7T + n - = 0.432 ± 0.007 and B^^^o^o = 0.217 ± 0.008. Using these values one gets for 
the ratio of charged to neutral decay widths, 

1 . -/ = 1.99 ±0.08 , 2.6) 

in perfect agreement with isospin symmetry conservation. 

3. Large- N c Chiral Perturbation Theory prediction 
3.1 Framework and lowest order calculation 

Large- Nq Chiral Perturbation Theory is the effective field theory of QCD in the chiral and 
large- Nc limits ||. In the large- Nc limit the U(1)a anomaly is absent and the pseudoscalar 
singlet 771 becomes the ninth Goldstone boson associated with the spontaneous breaking of 
U(3)l x U(3)r down to U(3)v- Both chiral and large- Nq corrections are treated pertur- 
batively. The effective Lagrangian is thus organized as a simultaneous expansion in powers 
of momenta, quark masses and l/Nc f| |32]| , 

C cS = £«» +£(*)+ £(2) + ... , (3.1) 
where the contributions of order 1, 5, <5 2 , ... follow the ordering of the series with 

8^ = 0(VS) , m q = 0(5) , l/Nc = 0{5) . (3.2) 
At lowest order, the = 0(e) ) has the form 3 ({A) stands for the trace of A) 

£(0) = il^Q^u) + fl {U i x + x tc/) _ l 2 M$ v l , (3.3) 

where x = ^BqM, Bq ~ 0(N C ) is related to the quark condensate, Ai = diag(m u , m^, m s ) 
is the quark-mass matrix, / ~ 0(\/Nc) is the pion decay constant in the chiral limit, Mq 
is the U(1)a anomaly contribution to the rji mass, and U = u 2 = exp (iV2$/f) with 



+ > + > - + K+ \ 



(3.4) 



In this case, the mathematical states rfe = (r/8,^1) are related to the physical states rfe 
(rj,rf) by 

(cos9p s'mOp 
— smOp cos Op 





(3.5) 



3 For the purpose of our work external sources are not required. The Lagrangian including them can be 
found in Ref. 101. 
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where Op is the 77-77' mixing angle in the octet-singlet basis at this order. At next-to-leading 
order, the part of = 0(5) relevant to the analysis of 77' — > r/mr is 

4'^- = L 2 (d^d v Ud^d v U) + (2L 2 + L 3 )(d^d»Ud u rfd»U) 

+L 5 (d fl U^d^U(U^ X + X ] U)) + L 8 (U^xU^x + x ] u x ] u) 



+-77-^771(^771 



(3.6) 



where the low-energy constants (LECs) L 2 , L3, L§, Lg are of 0(Nc), while Ai,A2 are of 
0(1/Nc). Ai and A2 only influence the singlet sector and can be attributed to OZI-rule 
violating contributions. Now, the relation between the mathematical and physical states 
cannot be simply expressed in terms of a mixing angle. The precise connection is found in 
Eq. (iBlTD . 

At lowest order (LO), i.e. using the in Eq. the amplitude is 0, |§ 0, 

l|, ID, |4] 

2 



A^ChPT I 



2\/2cos(20 P ) - sin(2(9p) 



6/ 2 



(3.7) 



where f = (the pion decay constant) at this order. For Op = (—13.3 ± 0.5)° |35[ and 
f n = 92.2 MeV, the branching ratio thus obtained is only 3% of the measured value. The 
main reason for this difference is the unexpected appearance of m 2 in the amplitude ( |3.7| ) 
that makes it to vanish in the chiral limit. Although in general the LO is able to provide 
a suitable approximation to the hadronic amplitude, e.g. the 77-77' mixing [17, |33| or the 
77(77') — > 77 decay [17, 34, 36, 37, I 



iq], in some few cases, like the anomalous magnetic 
moment g — 2 |3J| or the 77 — > 7r°7r° scattering [40], the LO contribution is absent (or very 
small) and the first (and dominant) contribution comes from higher orders. 

3.2 Next-to-leading order calculation 

At next-to- leading order (NLO), it is convenient to express the amplitude, from the point 

and 



of view of large- Nq, in terms of OZI-allowed (Ai 



and OZI-suppressed (M. 



r/ s riqiriT 



M 



T] s r] a 1TTT, 



contributions 4 , 



7]' — >riTTTV I NLO 



CqqM.T] q ri q -KTT 



+ C SS M 



(3.8) 



where the coefficients c qq , c sq and c ss (q stands for the u and d quarks and s for the s 
quark) are universal and encode the 77-77' mixing at next-to- leading order (see App. |B| for 
details) , 

f 



- qq 



3/ 2 / 2 cos 2 (# 8 -0 o ) 



2/| sin(20 8 ) - / 2 sin(2# ) - 2V2/ 8 /o cos(# 8 + ] 



4 The octet-singlet states (?78, 77i) are related to their quark-flavour counterparts {r) q ,rjs) by (778,771) = 
(77, - V^V s ,V2i] q + r/ s )/V3. 
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f- 



sq 



3/|/ 2 co S 2( 

f 2 



V2fl sin(20 8 ) + V2f 2 sin(20 o ) + / 8 / co S (9 8 + Q ) 



mfo cos 2 



/I sin(20 8 ) - 2/ 2 S in(20 o ) + 2a/2/ 8 /o cos(fl 8 + 0„) , (3.9) 



and 



VqVqKTr 



1 

7 1 



70; 



2(3L 2 + L 3 ) (s 2 + t 2 + y 2 



2L, 



f 2 

J 7T 



f' 2 

(to 2 , + to 2 + 2m 2 ) to, 2 + 



- m / 

24L 8 
f 2 



2m*) 



-m„ 



+ 3 A 2 < 



^2 
3/ 2 



A 2 m, 



2 

7T J 



A* 



ri s ri a -K-K 



(3.10) 



(3.11) 



where these amplitudes are independent of the way r] and rj do mix. Taking into account 
both the contributions of lowest and next-to-leading order, the amplitude for rf -4 r]irir is 
finally given by 

~m| 2(3L 2 + L 3 ) 
2 



ChPT 



-qq 



P 



f2 
J - 



(s 2 +t 2 + u 2 



2ml) 



2L5 
f 2 

J 7T 



(m 2 , + m 2 + 2m 2 ; 



m„ 



24L 8 



/ 2 * ' 3 



m 7r + o^2m n 



' *q V% . 2 



(3.12) 



As seen, the only contribution not proportional to m 2 and therefore dominant is the 3L 2 + 
L3 term. Indeed, this term is larger than the LO prediction which in turn is larger than 
the sum of the NLO L^,Lg and A 2 terms. Consequently, the branching ratio is expected 
to be of the form (0 + 6 m 2 + crn^) 2 ~ a 2 + 2 a bm 2 , where a corresponds to the 3L 2 + L 3 
term and b and c are the suppressed contributions. 

For the numerical analysis, we use fg = 1.28/^, fo — 1.25/„-, # 8 — —20°, 9q ~ —4° and 
A 2 ~ 0.3 from Ref. EL and = 137.3 MeV, m v = 547.9 MeV and raw = 957.8 MeV 



from Ref. |23|]. For the LECs in the large- Nc limit we find in the literature three different 
sets of values gathered in Table ^. All of them come from resonance-exchange estimations 
with only one resonance per channel [42, 43, 44 1. In the simplest scheme one obtains the set 
of LECs identified as setl f42fl . The values of the LECs are slightly modified after including 
QCD-inspired assumptions of high-energy behaviour 5 . When an unsubtracted dispersion 
relation for the pion vector form factor is used one obtains set2 [H] . When constraints on 
the scalar form factor are also used one gets set3 El. 

The predicted branching ratio for the different sets of LECs is found to be 14% (setl), 
9% (set2) and 64% (set3). As seen, the next-to-leading results are in general one order 



5 Using a model with a finite number of resonances instead of an infinite one, which would correspond 
to the full large- Nc theory, may introduce problems in the short-distance matching and large discrepancies 
in the prediction of LECs fia| . 
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setl fh| 


setg | 


43 


1 


set3 \ 


44 


1 


L 2 • 10 3 


1.2 




li 






li 




L 3 • 10 3 


-3.0 




4.9 




4.3 


(3L 2 + L 3 ) • 10 3 


0.6 




0.5 




1.1 


L 5 • 10 3 


1.4 




1.4 




2.1 


L 8 • 10 3 


0.9 




Oi 


) 




Oi 





Table 2: Values for the three different sets of low-energy constants discussed in the text. 



of magnitude bigger than the lowest order prediction after Eq. (3.7), thus indicating that 
our approach for the description of rf — > r/TTir moves into the right direction. However, 
due to the broad range of values obtained for the predicted branching ratio, caused by the 
numerical difference in the 3L2 + L3 combination, it seems to be more appropriate to fix the 
value of this dominant contribution from the measured branching ratio while keeping the 
remaining LECs to the values shown in Table |2| for the corresponding set. Since 3L 2 + L3, 
L5 and L% are similar in setl and set2, and different with respect to set3, we will use for 
the purpose of comparison only setl (without QCD-inspired constraints) and set3 (with 
QCD-inspired constraints). In Table ||, the outcome for the 3L2 + L3 combination obtained 
in this way are displayed. The values seem to be in line with set3. 

Besides the branching ratio, the study of the Dalitz plot distribution provides relevant 
information. The squared amplitude \M\ 2 is expanded in terms of the variables X and Y 
(see Sec. ^ for details), where Y is a linear function of s, the tttt invariant mass, and X 
appears always in the form of cos^ = Xf(Y), with 6 n the angle of p„.+ with respect to p„ 
in the tttt rest frame. This explains the reason for the absence of odd powers of X under 
the assumption of C-invariance. Interesting facts concerning the Dalitz plot parameters are 
the following. First, the amplitude ( p. 12 ) in the chiral limit (with m^m^i 7^ 0) becomes 

lim^M^X^Tr = -^ ——7T^ — m v/ (m v > - m v ) 2 (6m v - m v >X 2 ) . (3.13) 



111 



This means that when A4 is expanded in X and Y around the center of the Dalitz plot the 
terms with powers of Y will be proportional to m n and because of that suppressed. As a 
result, we find that the terms in the expansion of order Y are similar in magnitude to those 
of order X 2 . Then, one would expect for the Dalitz plot parameters the following hierarchy: 
a ~ d 3> b ~ «2i ~ ^40- For this reason, we propose to extend the standard Dalitz plot 
parameterization used by experimental analyses (based on a, b and d alone) and include 
the additional K21 and K40 parameters for consistency. Second, certain combinations of 
parameters turn out to be independent of the chiral couplings, being just functions of the 
pseudoscalar masses. In particular, we obtain 

a _ «2i _ §rn n {m v - m 7r )(m^ / - (m v + 2m 7r )(3m r? / - m v - m n )) _ ^ ^ 
d 2^40 m r] i(jn ri + 2m ir ) 2 (m r) ' — — 2m- K ) 
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Our predictions for the Dalitz plot parameters are summarized in Table ||, together with 
the reported values from GAMS-4-7T || and VES || based on analyses of 7/ — > ?77T 7r and 
7/ — > 7]tt + tt~, respectively. The parameters are almost identical for both sets thanks to 
the dominance of the 3L2 + combination over the tiny contributions from L5, and 
A2 terms. When compared to the measured values the d parameter is in the correct range 
whereas the a and b are clearly incompatible. This points out the need for including higher 
order effects. Final state interactions originated from the rescattering of the two pions are 
the main source of such effects. 



Parameter 


setl 


set3 


GAMS-4tt | 


I 


VES £ 


I 


a[Y] 

b[Y 2 } 

d[X 2 } 

K2l[X 2 Y] 
«4o[A 4 ] 


-0.303 
+0.001 
-0.089 
+0.014 
+0.002 


-0.284 
-0.001 
-0.084 
+0.012 
+0.002 


-0.066 + 0.016 
-0.063 ± 0.028 
+0.018 ± 0.078 


-0.127 + 0.018 
-0.106 + 0.031 
-0.082 + 0.019 


(3L 2 + L 3 ) ■ 10 3 


1.0 


0.9 







Table 3: Dalitz plot parameters for setl and setS and the measured values (errors are added in 
quadrature) from rj' r]ir ir (GAMS-47r) and rf — > 7771-+ 7r~ (VES). The combination ZL2 + £3 is 
extracted for each set from the experimental branching ratio. 



3.3 7T7T final state interactions 

As stated, a further improvement in the large- Nc approach would be to calculate the next- 
to-next-to-leading (NNLO) effects. In this case, loops come into play and the LECs are 
scale-dependent. A detailed calculation of the loop effects would account for the related 
final state interactions. In view of the exploratory nature of our analysis this calculation is 
beyond the scope of the present work. However, we can estimate these rescattering effects 
by means of a unitarization procedure. Different from scattering processes, in 7/ — > 777777 all 
the channels contribute to the unitarity relation (see App. [A.2| for details). Nonetheless, 
the contributions from rjir rescattering effects {t- and u-channel) are negligible. This was 
also remarked in a previous non-relativistic effective field theory study pTj], where they 
found a very small value for the tjtt scattering length. Accordingly, we only consider s- 
channel unitarity due to 7T7T final state interactions. This unitarity constraint is better 
expressed in terms of partial waves (see Eq. ( |A.9| )). The constraint incorporates the 7T7t 
1 = partial- wave amplitudes Tj, which for J = 0, 2 are found in App. A.l. 



For r( —7- 777T7T the relevant partial waves are J = 0, 2 (higher partial waves are sup- 
pressed by phase space) since I(tttt) = and I + J = 0, 2 . . . must be fulfilled. Hence, we 
compute the S- and D-wave projections of the amplitude ( 3.12 ) as specified by Eq. ( |A.2j ). 
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The result is 

1 



M (s) 



3m 2 + 4 ^ 3 ^ 2 J" ((s - m^/)(5s + (m 2 , - 8m 2 ) + 2m 2 im 2 /s) 

/•7T 



192vr \ p 

—Asm 2 + m^(2m^ — m 2 + 10m 2 ) + 2m^m^(2m^ — m 2 )/s) 
-12ml (~^( m v' + m ^ + 2 m l) - m l \ + 4A 2 m 2 



+£|2V2A 2 m 2 j , (3.15) 

M2{S) = 2^ C f 3L2 tf L3 {s ~ 4m ' } (S " 2(m '' + + {m "' ~ m ' )2/s) • (3 - 16) 

The predicted branching ratio is completely dominated by the S- wave (scalar component), 
while the Dalitz plot parameters are not. The reason is that the variable X always appears 
through cos 6^ and, as a consequence, if only the 5-wave is considered the terms propor- 
tional to X are absent and cf[AT 2 ] = /t 2 i[AT 2 y] = K4o[X 4 ] = 0. The first contribution to 
the former three parameters comes from the D-wave, which turns out to be essential for 
their precise determination. For the a and b parameters the contribution of the D-wave is 
clearly less important than the scalar component. 

Once the partial wave amplitudes are known, one can proceed to unitarize them. 
We use two different procedures (see App. |A.2j for details), one based on the if -matrix 
formalism in Eq. ( |A.10[ ) and the other inspired by the NjD method in Eq. ( [A.12 ). The 
final rf — > rjirir unitarized invariant amplitudes are 

M(s,t,u)\ K . matlix = V 32vr(2J + l)P J (cos 6 W ) - ^/j^™ ' ( 3 - 17 ) 
in the simplest if -matrix formalism, and 

M(s,t,u)\ N/D = ^M2J + l)Pj(co S 9 n ) , (3.18) 

in the N/D method. 

Once these amplitudes are computed, the corresponding Dalitz plot parameters can 
be extracted. The results within the if -matrix unitarization procedure are shown in Table 
^ The combination 3L 2 + L3 is again fixed for each set of LECs from the experimental 
branching ratio. Notice that when compared with the results in Table [3L the parameters 
only carrying powers of X, that is d and K40, do not change, thus indicating that the 
unitarized version of the D-wave does not produce any quantitative modification. The 
reason is that the 7j 7T7r-scattering amplitude entering into the unitarization of the D-wave 
emerges at next-to-leading order in the chiral expansion. On the contrary, the parameters 
carrying powers of Y are substantially modified because of the unitarization of the 5-wave 
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Parameter 


setl 


set3 


GAMS-4vr | 


! 


VES % 


I 


a[Y] 
b[Y 2 ] 
d[X 2 ] 
K2i[X 2 Y] 


-0.267 
-0.018 
-0.089 
+0.010 
+0.002 


-0.248 
-0.020 
-0.084 
+0.009 
+0.002 


-0.066 + 0.016 
-0.063 ± 0.028 
+0.018 ± 0.078 


-0.127 + 0.018 
-0.106 + 0.031 
-0.082 + 0.019 


(3L 2 + L 3 ) • 10 3 


1.0 


0.9 







Table 4: Dalitz plot parameters obtained from the amplitude in Eq. (3.17) using the JT-matrix 
formalism. The experimental values are also shown for comparison. 



(the main contribution to Tq is in this case the leading order). When compared to the 
measured values, the parameters a[Y] and b[Y 2 ] have changed in the right direction, but 
not enough, and d[X 2 ] still agrees. 

For the unitarization inspired by the N/ D method (a sophisticated version of the K- 
matrix formalism), the constant C involved in the Bq(s, m 2 , m 2 ) integral of Eq. ( |3.18 ), see 



also Eq. (A. 11), must be fixed from experiment together with the combination 3L 2 + L3. 
In order to fix C we use the better measured Dalitz plot parameter, a[Y], which is either 
-0.066(16) [§ or -0.127(18) f§. This means that the parameter ranges -0.145 < a < 
—0.050, or equivalently a = —0.098(48), in accordance with both experiments. We keep 
using the experimental branching ratio to fix 3L 2 + L3 . The results are shown in Table [5[ 
The parameters b and d are now in line with the measured values. 



Parameter 


setl 


set3 


GAMS-4tt § 


VES § 


b[Y 2 } 

k 21 [X 2 Y] 
Kw[X 4 ] 


-0.051(1) 
-0.101(8) 
+0.004(2) 
+0.003(1) 


-0.050(1) 
-0.092(8) 
+0.003(2) 
+0.002(1) 


-0.063+ 0.028 
+0.018+0.078 


-0.106+ 0.031 
-0.082+ 0.019 


(3L 2 + L3) ■ 10 3 
C-constant 


1.1(1) 
-2.0 < C < 0.4 


1.0(1) 
-1.4 < C < 0.9 







Table 5: Dalitz plot parameters obtained from the amplitude in Eq. ( 3.18 ) inspired by the N/D 
method. The range of the constant C is fixed from a[Y] = —0.098(48). The errors account for the 
propagation of the uncertainty in a[Y\. 



To sum up, rescattering effects are seen to be important in r( — > -qinr. The unitarization 
procedure give us an idea of what is missing when the amplitude is computed only up to 
NLO in the large- Nc ChPT expansion. From the results in Table || to those in Table ||, there 
is a clear improvement on the predictions for the Dalitz plot parameters as compared to the 
measured values. Finally, the different parameters agree well within experimental errors. 
The next logical refinement would be to incorporate the lowest lying set of resonances and 
check wether they are the sole responsible for such improvement. We postpone this analysis 
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to the next section. 

For the benefit of future experimental analyses, we display in Fig. |], the Dalitz plot 
distribution of rf —> rjirn using Eq. ( p. 12 ), and in Fig. ^J, the corresponding invariant mass 
spectra from Eq. (|3.12| ) — blue dashed line — and Eq. ( |3.18| ) — solid orange band. In both 
figures, we use set3 for the LECs. 



> 

01 

O 




U 




-1.0 -0.5 0.0 0.5 1.0 



Mj (GeV 2 ) 



Figure 1: Dalitz plot distribution of rf —> rj-K-K using Eq. (3.12) supplemented by rescattering 
effects through Eq. ( |3.1Sj ), in terms of the invariant masses M^ ff and MjL (left) and the kinematical 
variables X and Y (right). Larger values are shown lighter. 
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0.08 0.10 0.12 0.14 0.16 



Figure 2: M™ (left) and (right) invariant mass spectra for the differential branching ratio. 



The tree-level large- Nc ChPT prediction from Eq. (3.12) — blue dashed line — is compared to its 



unitarized counterpart via Eq. (3.18) — solid orange band 



4. Resonance Chiral Theory prediction 

As we discussed in the previous section, large- Nq ChPT at NLO provides a successful pre- 
diction for the Dalitz parameters. However, the comparison of Tables || and |5] shows that 
final state interactions effects (NNLO in the large- Nc ChPT counting) are not negligible 
and must be properly incorporated. One may wonder then whether other possible NNLO 
ChPT contributions could be relevant, such as higher dimension local terms in the chiral 
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Lagrangian (e.g. 0(p e ) operators). Previous studies have found that these local contri- 



butions can be saturated by tree-level resonance exchanges j|, [46], |47[], where the heavy 
meson propagators resum higher chiral orders. Thus, we will make a thorough estimate of 
the impact of these higher chiral orders by observing the contribution from intermediate 
resonant states. 

4.1 Framework and lowest order calculation 

Resonance Chiral Theory (RChT) is a description of the Goldstone and resonance inter- 
actions in a chiral invariant framework [Q, fl(], ^7|]. The pseudo-Goldstones enter through 
the exponential realization u = exp (i^/v^/)- As the standard effective field theory mo- 
mentum expansion is not valid in the presence of heavy resonance states, RChT takes 
the formal l/Nc expansion as a guiding principle |48|]. Thus, at large- Nc, the mesons 
get classified in U(nj) multiplets, with nj = 3 the number of light-quark flavours. The 
pseudo-Goldstone multiplet gathers then the octet of light pseudo-scalars (n,K,r]s) and 
the chiral singlet rji, conforming a U(3) nonet <£. For convenience, the Lagrangian can be 
then organized according to the number of resonance fields in the interaction terms, 

C RxT = C GB + C Rl + C R ^ + £ R i R j R >° + ■■■ , (4.1) 

where Ri stands for the resonance multiplets and the first term in the r.h.s. of the equation 
contains the operators without resonance fields, 

£ GB = ^( U / + X+ ) , (4.2) 

with = iv)d^Uv) = ujt (external sources are not included) and x+ = u^X u ^ + u\^u. For 
scalar resonances, 5'(0 ++ ), the second term in Eq. ( f4.1|) is given by the operators jl(| 

C s = c d (Su^) + c m (S X+ ) , (4.3) 

where only the terms with a chiral tensor of order p 2 have been shown. The resonance 
Lagrangian will be supplemented with kinetic and mass terms for the resonance fields. 
The Lagrangian (4.1) also contains operators with vector meson fields, but these do not 



contribute to rf — > ijtttt because of G-parity conservation. 

A priori, one should also consider resonance operators in £r, x t with chiral tensors of 
arbitrary order 6 . For the case of scalar resonance operators, one can prove through meson 
field redefinitions that the most general scalar-pseudoscalar-pseudoscalar interaction in the 
chiral limit is provided by the c<j term in Eq. (|4.3j) |49|. On the other hand, although nothing 



6 The infinite tower of mesons contained in large- Nc QCD is often truncated to the lowest states in each 
channel with successful predictions for the 0(p 4 ) and 0(p 6 ) low-energy constants. However, it has been 
pointed out in Ref. E3| that large discrepancies may occur between the values of the masses and couplings 
of the full large- Nc theory and those from descriptions with a finite number of resonances. 
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prevents the presence of £ GB operators of order higher than p 2 , within the large- Nc limit, 
a bounded behaviour of the rf — > rjirn partial-wave amplitude at high energies requires the 
absence of 0(p 4 ) operators and higher in £ GB . Consequently, the sum C GB + C s contains 
all the pieces of the RChT Lagrangian that one would need to describe the rf — > rjtnr decay 
in the large- Nc limit. Needless to say, the U{1)a anomaly contribution to the singlet mass 
must be also taken into account by means of the Lagrangian [47] 



C 



P 



m 



Ml In 2 (det u) 



1 



(4.4) 



3 u v ' 2 

This contribution, which is not fixed by symmetry requirements alone and depends crucially 
on the dynamics of instantons, is formally next-to-leading order in the large- Nq expansion. 
However, it is essential for the correct description of the rj and its decays. Apart from this 
particularly large effect, any other NLO effect in 1/Nc is expected to be suppressed. 

Concerning the 7/ — > r]7nr calculation in the framework of RChT, we follow Sec. 3.2 and 
express the amplitude in terms of OZI-allowed and OZI-suppressed contributions. At lowest 
order in large- N c , the latter vanish, = M Vsns irK = 0, and Mrf^rpnt = c qq Mn qriq n-K 

is given by 



M 



RChT 

77'— trjTTTT 



^qq 
P 



2 



+ 



1 (c<i(s 



mz 



ml) + 2c m ml){c d {s - 2m\) + 2c m ml] 



f- 

J 7T 



Ml 



1 (c d {t 



m 



ml 



+ 2c 

in 



mz 



ml) + 2c m m 2 ) 



+ 



J 7T 

1 {c d (u-m 2 , 



M 2 



+ 2c m ml)(c d (u 



m~ 



+ 2c m ml) 



J 7T 



Ml 



u 



(4.5) 



where one has the contribution of the 1 = 1 scalar resonance ao in the t- and n-channel 
and the 1 = scalars a and fo in the s-channel. Ms is the mass of the scalar multiplet. 



Although the Lagrangian (O) does not include scalar mass splitting nor a-fo mixing, these 
can be incorporated in the theory in a straightforward way through an operator of the form 



|53]. Without loss of generality, one may consider different masses for the 



resonances in the scalar multiplet and a mixing scheme for the a and /q. In the heavy 



scalar mass limit, that is p 2 ~ m 2 p ^ Mj, the amplitude in Eq. 



M 



RChT 
r} 1 — >r]7nr 



■qq 



P 



m^ 
2 



+ 



PM 2 



[s 2 + t 2 + u 2 



m r 



■jp^farf +m v + 2m n )m n + 



m r 



m„ 



becomes 



2m! 



(4.6) 



which turns out to be analogous to the large-iVc ChPT amplitude in Eq. fljTI^ up to 
subleading contributions in 1/Nc- After using the low-energy constant relations 3L2+-L3 = 
c 2 d /2Ml, L 5 = c d c m /Ml and L 8 = c 2 m /2Ml @, and the large-iVc result Ai j2 = 0, the two 
amplitudes are seen to be formally identical. 
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The largest contribution in Eq. fll.5| ) to the branching ratio comes from the q terms, 
which are proportional to the external momenta. Everything else is proportional to m 2 , 
being suppressed. Nevertheless, even though the decay rate is essentially proportional to 
c d , the interference terms in the squared amplitude are not numerically negligible and must 
be kept at the current level of experimental precision. In addition to the branching ratio, 
the study of the shape of the differential decay width provides important information. For 
instance, the amplitude in Eq. ( |4.5[ ) only depends on X 2 in the chiral limit, 

lim uRChT _ c qq 6c 2 q + a 2 X 2 

where a n and (3 n are some polynomial combinations of just and Ms- This pure 

X 2 dependence was also found before in the large- Nc ChPT framework but with simpler 
polynomial combinations of pseudo-Goldstone masses. 

For the numerical analysis, we use the same input data we utilized in Sec. 3^2 for the 
mixing parameters and pseudoscalar masses, together with Ms = 980 MeV for the scalar 
multiplet mass. For the scalar couplings Cd and c m , the latter, which is less relevant since 
it always appears through chiral suppressed contributions, will be fixed via the high-energy 
scalar form factor constraint 4crfC m = f 2 (indeed, 4crfC m ~ f 2 is taken) j54j. Different 
values for these couplings are reported in the literature. From the ao(980) decay width, 
Cd = 26±9 MeV and c m = 80±21 MeV 152]; from low-energy constant resonance saturation, 



c d = 30 ± 10 MeV and c m = 43 ± 14 MeV [||, ||]; from the I =1/2 Kir s-wave scattering 
amplitude analysis, Cd = 25 MeV and c m = 77 MeV [55], and from the / = 3/2 Kit s-wave, 
Cd = 13 MeV and c m = 85 MeV [55|. Notice that all these results obey rather well the 
theoretical constraint 4crfC m = f 2 . 

We mentioned that the largest contribution to the branching ratio comes from the Cd 
terms. Indeed, it is to a large extent proportional to the fourth power of Cd- Therefore, 
owing to the large uncertainty in q, we prefer to fix the value of this coupling from the 
experimental branching ratio and predict the Dalitz plot parameters. In doing so, one gets 
Cd — 28 MeV, in agreement with most of the reported values. A first prediction of the 
Dalitz plot parameters is shown in the second column of Table ^. While a and d are in 
nice agreement with the VES results, b is very far. We also calculate certain combinations 
of these parameters that were found to be independent of the chiral couplings in the large- 
Nc ChPT framework. In the case of RChT, even though the cancelation of couplings 
is not complete, there is not a sizable dependence on the value of Cd- In particular, one 
finds a/d = 2.1 and K21/2K40 = 2.4, to be compared with a/d = ^21/2^40 = 3.4 in 
Eq. ( |3.14| ). The RChT result for the first ratio agrees better with the measurement from 
VES, a/d| exp = 1.55 ± 0.42, than the large-iVc ChPT prediction. 

In addition to the contribution of scalar mesons, we also consider the impact of J = 2 
resonances. In the previous large- Nq ChPT analysis, the J = 2 partial wave was shown to 
be crucial to properly recover the Dalitz plot parameters. Still, as the mass of the lightest 
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Parameter 


J = 


J = 0,2 


GAMS-4vr § 


VES § 


a[Y] 

b[Y*] 

d[X 2 ] 

K2l[X 2 Y] 
«40[^ 4 ] 


-0.119 
+0.001 
-0.056 
-0.006 
-0.001 


-0.160 
-0.001 
-0.070 
-0.003 
-0.001 


-0.066 ±0.016 
-0.063 ± 0.028 
+0.018 ±0.078 


-0.127 ±0.018 
-0.106 ±0.031 
-0.082 ±0.019 


c d (MeV) 


28 


26 







Table 6: Dalitz plot parameters obtained from the RChT amplitude in Eq. (4.5) (second column) 



and with the addition of J = 2 resonances via Eq. (4.8) (third column) 



tensor multiplet is roughly Mf 2 = 1.2 GeV, one may just consider its leading effect in 
1/Mj 2 rather than the whole non-local resonance propagator structure. Thus, it induces a 
contribution that is identical to the 3i>2 + £3 term in Eq. ( |3.12| ), 



c qq 2(3Ll + LT)^ 2 + t2 + u2 



in, 



2<) 



P ft 

The J = 2 resonance contributions to the 0(p 4 ) LECs, 3L\ + Lj = gjj3M 



(4.8) 



0.16-10- 



were estimated in Ref. [56], after imposing high-energy constraints on 7T7r-scattering. When 
the tensor contribution is taken into account, the value of the scalar coupling obtained from 
the experimental branching ratio is ~ 26 MeV. The predicted Dalitz plot parameters 
in this case are shown in the third column of Table ^. The values for a and d are still in 
reasonable agreement with VES, but b continues to be very far. 



4.2 Resonance width effects 

When comparing Tables ^ and [6|, one can see that RChT at tree-level improves considerably 
the prediction for the a parameter from NLO large- Nc ChPT. Nonetheless, in both cases, 
one finds a far too small value for b. This issue was solved in large- Nq ChPT by considering 
the FSI impact from meson loops through a convenient unitarization scheme. It is then 
desirable to incorporate loops in the resonance description in order to refine our prediction, 
the same we did in the ChPT framework. The first obvious FSI effect is that the scalar 
mesons gain a non-zero width. In particular, the scalar isoscalar a meson becomes broad. 
Although a detailed study of these NLO effects in 1/Nc is beyond the scope of this work, 
we perform some estimates of the theoretical uncertainty associated with this large- Nc 
description. 

A first estimate is provided by the inclusion of a cr-/o(980) splitting and a self-energy 



in the s-channel propagator of the a, by means of the substitution in Eq. (4.5), 



Ml 



sin 2 4>s 

M 2 f -s 
to 



+ 



cos 4>s 



Ml 



s - c a s k B (s,ml,ml) 



(4.9) 
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with Mf = 980 MeV, 4>s = —8° Q and the parameters M a and tuned such that one re- 
covers the right position for the a pole, M| ole = 441+^ 6 MeV and lT le = 544^5 MeV @. 
The function Bq is the subtracted two-point Feynman integral written in Eq. ( A, lip with 
C = 2 and p(s) = yl — 4m^/s. The power behaviour = produces an unphysical bound 
state in the first Riemann sheet very close to the tttt threshold from below, which enhances 
unnaturally the amplitude and leads to ~ 10 MeV. This case seems to be disfavoured 
from the phenomenological point of view and is discarded. For k = 1, the amplitude repro- 
duces the a pole and fixes Cd — 26 MeV. Power behaviours equal or higher than k = 2 are 
unable to generate the a pole at the right position, and are also discarded. The predicted 
Dalitz plot parameters for k = 1 are shown in Table |7[ As seen, the value of b moves into 
the right direction, a increases with respect to its value without including cr-width effects, 
and the rest of parameters present small variations (see Table ^ for comparison) . One can 
take also into account the /o(980) and ao(980) widths in a similar way 7 , but their effects are 
seen to be clearly subdominant as compared to the impact of a broad a and they slightly 
change the Dalitz plot parameters obtained before. 



Parameter 


J = 
+o width 


J = 
+ (cr, fo,a ) widths 


GAMS-4^ § 


VES § 


a[Y] 


-0.161 


-0.168 


-0.066 ± 0.016 


-0.127 ±0.018 


b[Y 2 } 


-0.034 


-0.035 


-0.063 ± 0.028 


-0.106 ±0.031 


d[X 2 } 


-0.055 


-0.056 


+0.018 ± 0.078 


-0.082 ±0.019 


K 2 i[X 2 Y] 


-0.005 


-0.006 








-0.001 


-0.001 






c d (MeV) 


28 


26 







Table 7: Dalitz plot parameters obtained from the RChT amplitude in Eq. (4J3) supplemented by 
cr-width effects (second column) and with the addition of fa- and ao-width effects (third column). 



The estimate of the FSI effects provided here is obviously model dependent, as we have 
introduced ad hoc a splitting and a self-energy for the scalar multiplet. The splitting can 
be easily introduced through the corresponding terms in the Lagrangian, an issue studied 
in detail in Ref. |53| . On the other hand, the resummation of the one- loop self-energy 
effects is also justified, even for the case of the broad a. Higher order effects as multimeson 
channels are negligible below 1 GeV and the one-loop amplitude seems to provide the 
relevant information in the physical range under study. In any case, a self-energy must be 
appropriately resummed in the neighbourhood of the resonance pole f60| , 61]. Nevertheless, 



as we are aware of the level of model dependence of the present procedure, in the next 



7 For the fo and ao scalar mesons, the important piece of the self-energy is its imaginary part, being the 
real part of its corresponding logarithm almost negligible in comparison with the dominant contribution 
Mg — k 2 (with k 2 — s, t, u, depending on the channel). In any case, we consider the full self-energy function 



as done in Eq. (4.9). The function B (k 



i 2 ) appearing in the t- and w-channel can be found in Ref. 



- 18 - 



subsection we use the N/D unitarization method to incorporate meson rescattering effects 
in a less dependent way and compare both approaches. 



4.3 7T7T final state interactions 



As in large- Nq ChPT, here we consider FSI by demanding unitarity in the tttt channel. We 
have already seen that this is by far the most important rescattering effect due to the little 
impact of the ao width (see Table [?]) . We use again the N/ D method to unitarize the tree 
level RChT amplitude in Eq. ( |4.5| ) trough the expression in Eq. ( grp . In this manner, 
provided the absorptive cuts in the t- and u-channel are neglected, only the 7T7r-scattering 
amplitude in the s-channel has to be specified and not the precise way the resonances get 
their width. The scalar coupling a and the constant C appearing in the Bq integral of 
Eq. ( |A.11| ) will be fixed from the experimental branching ratio and the a parameter. The 
predictions for the other Dalitz plot parameters are given in Table |8[ We consider these 
results as the most reliable RChT estimate. The comparison of the third column in Table ^ 
and the second column in Table || shows a fair agreement between the two FSI approaches, 
resonance propagator modification and N/D unitarization. The largest difference is in the 
a[Y] parameter, which is taken as an input in the N/D method. 



Parameter 


J = 
+N/D unitar. 


J = 0,2 
+N/D unitar. 


GAMS-4-7T § 


VES § 


b[Y 2 ] 
d[X 2 ] 

K2l[X 2 Y] 
K40[X 4 } 


-0.030(1) 
-0.057(1) 
-0.010(1) 
+0.001(1) 


-0.033(1) 
-0.072(1) 
-0.009(2) 
+0.001(1) 


-0.063 ± 0.028 
+0.018 ± 0.078 


-0.106 + 0.031 
-0.082 + 0.019 


c d (MeV) 
C-constant 


27(1) 
2.1 < C < 3.4 


25(1) 
1.4 < C< 2.7 







Table 8: Dalitz plot parameters obtained from the RChT amplitude in Eq. (4.5) supplemented by 



rescattering effects via Eq. (3.18) (second column) and with the addition of J = 2 resonances (third 
column) . 



To summarize, if one remains at tree level there is no way to increase the prediction 
for the b[Y 2 ] parameter from a value very close to zero to the size of measured values. 
The same was found in large- Nc ChPT. This situation is greatly improved as soon as one 
incorporates the rescattering effects, see Table ||. Now, the b parameter is compatible with 
the GAMS measurement and not far from the VES result. On the other hand, the d[X 2 ] 
parameter seems to be rather stable, both in large-Nc ChPT and RChT, either at tree 
level or including rescattering effects. It remains compatible with the various experimental 
values, in particular, with the most precise measurement reported by VES. 

Finally, we display in Fig. |3|, the Dalitz plot distribution of if — > r]7nr using the RChT 
amplitude in Eq. (JOJ) supplemented by rescattering effects through Eq. (3.18), and in 
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Fig. ||, the corresponding invariant mass spectra from Eq. (4.5) — blue dashed line — and 
Eq. ( gig) —solid orange band. In Fig. ||, we use as input the values shown the third 
columns of Tables ^ and || for the lowest order calculation and its N/D unitarized coun- 
terpart, respectively. 



> 

01 

O 




o 

-1 



Mj (GeV 2 ) 




Figure 3: Dalitz plot distribution of rj — > ijTnr using Eq. (4.5) supplemented by rescattering 
effects through Eq. fl3.1S| ), in terms of the invariant masses and MjL (left) and the kinematical 
variables X and Y (right). 
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Figure 4: Mz^ (left) and M% n (right) invariant mass spectra for the differential branching ratio. 



The tree-level RChT prediction from Eq. (4.5) — blue dashed line — is compared to its unitarized 



counterpart via Eq. (3.18) — solid orange band 



5. Discussion and conclusions 

In this work, we have studied the hadronic decays rj' — > t/tttt in the frameworks of large- Nq 
Chiral Perturbation Theory, at lowest and next-to- leading orders, and Resonance Chiral 
Theory in the leading 1/Nc approximation. In both cases, we have also considered higher 
order effects such as 7T7t final state interactions through detailed unitarization procedures. 
The inclusion of finite-width effects in the case of RChT is also discussed. In addition, and 
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for the benefit of present and future experimental analyses, we have computed the Dalitz 
plot distribution and the differential branching ratio. 

In large- Nc ChPT, the next-to-leading order calculation in Eq. ( 0.12 ) represents a clear 
improvement with respect to the current algebra expression in Eq. (|3.7|) , leading order in 
large- Nc ChPT, and yields a prediction for the branching ratio of the right size. The RChT 



computation in Eq. (4.5) also constitutes an improvement over the current algebra result. 
However, the two approaches fail to reproduce the Dalitz plot parameters a[Y] and fo[5^ 2 ], 
both associated with s-channel exchanges and then related to isoscalar scalar resonances. 
This points out the importance of a meson effects in this channel. On the contrary, the 
contributions from the t- and tt-channel are less problematic. The predictions for the 
parameter <i[X 2 ] in the different frameworks are relatively consistent among themselves 
and agree with the experimental measurements. 



Parameter 


Large- N c ChPT 
+N/D unitar. 


RChT (J = 0, 2) 
+N/D unitar. 


GAMS-4tt § 


VES § 


a[Y] 


-0.098(48) + 


-0.098(48)+ 


-0.066 + 0.016 


-0.127+0.018 


b[Y 2 } 


-0.050(1) 


-0.033(1) 


-0.063 + 0.028 


-0.106+0.031 


d[X 2 } 


-0.092(8) 


-0.072(1) 


+0.018 + 0.078 


-0.082+ 0.019 


n 2l [X 2 Y] 


+0.003(2) 


-0.009(2) 






K i0 [X 4 } 


+0.002(1) 


+0.001(1) 






(3L 2 + L 3 ) ■ 10 3 
c d (MeV) 
C-constant 


1.0(1) 
-1.4 < C < 0.9 


25(1) 
1.4 < C < 2.7 







Table 9: Final results for the Dalitz plot parameters obtained in the frameworks of large- Nc ChPT 
(second column) and RChT with the addition of J = 2 resonances (third column), respectively, 
after the unitarization of the corresponding amplitudes through the N/ D method. The rf — ¥ rjmr 
experimental branching ratio and the Dalitz plot parameter a have been used as inputs^ to fix the 
relevant couplings and the constant of the unitarization formula. 



Our results can be further improved through an appropriate unitarization of the am- 
plitudes in order to account for the final state interactions in the tttt system. The best 
estimates are achieved with the N/D unitarization method. The final results for the Dalitz 
plot parameters using this method are summarized in Table [9|. We show the results for set3 
of low-energy constants in the case of large- Nc ChPT and including the effects of scalar 
and tensor resonances for the case of RChT. For the sake of comparison, we also present 
the measured values. All Dalitz plot parameters lie now within the experimental range. 
The differences between the predictions of the two frameworks in Table p| should be taken 
as the systematic uncertainty of our analysis. In any case, it is much smaller than the dis- 
persion exhibited by the measurements. Preliminary results from BES-III seem to provide 
much more precise determinations, a = -0.047 + 0.011 ± 0.003, b = -0.069 ± 0.019 ± 0.009 
and d = —0.073 ± 0.012 ± 0.003 [59). If confirmed, this would favour our predictions with 
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respect to other theoretical studies which predict, for instance, a very small and positive d 
parameter (see Table [l] and Ref. [^] ) . 

In conclusion, the large- Nq ChPT framework seems to be suitable for describing rj — > 
r/inr decays even though higher next-to-next-to- leading chiral orders turn out to be relevant. 
The RChT approach tries to cure this problem since part of the higher order local terms 
in ChPT are now resummed through the exchange of heavy resonances. We have seen 
that the ao(980) contribution is dominant but, at the same time, the a contribution is 
also essential to agree with experiment. However, we find that heavier states, which are 
not present in the original RChT framework, are not negligible. In particular, the lowest 
lying tensor resonances produce a noticeable effect in the RChT estimates. Therefore, 
although RChT resums higher chiral orders, at this level of precision, it requires a much 
more detailed knowledge of the resonance content and properties. On the other hand, 
in the large-iVc ChPT approach, the effect of any possible heavy state is encoded in the 
low-energy constants and hence the resonance spectrum does not need to be known in 
detail. For these reasons, we consider both frameworks, large-Ac ChPT and RChT, to 
be complementary, where each one has its own advantages and disadvantages. We hope 
future experiments to be able to distinguish the most convenient framework for the study 
of this and other rj decays. 
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A. Partial waves and unitarity 



A.l Partial- wave projection 

The projection of a given rf — > tjtttt invariant amplitude, Ai(s,t,u), into partial waves is 

M(s,t,u) =^32vr(2J + l)Pj(cos^)7Wj(s) , (A.l) 
J 

where Pj is the Jth Legendre polynomial, 9 W is the angle of p„.+ with respect to p^ in the 
7T7T rest frame, and 

^> = m A(.», m? X,.)^(,<,<)W , *^(~^(..*.«) , (A.2) 



with 

s(m 2 , + m 2 + 2ml - s - 2i) 
X(s,m 2 ,m 2 ,) 1 / 2 X(s,m 2 ,m 2 ) 1 / 2 ' 



cos « = ~ TT \ 2TT79T7^ 9 , > ( A - 3 ) 



^max(min) 2 



A(s,m2,m2,)V2 A ( S)m 2 ) m 2)l/2 

m v , +m v + 2m n — s± ■ 



, (A.4) 



and A(x, y, z) = x 2 + y 2 + z 2 — 2xy — 2xz — 2yz. For a 7T7r-scattering amplitude of isospin 
/, its decomposition into partial waves is 

T J ( S , t, u) = £ 32^(2 J + l)Pj(l + J^2)Tj{s) , (A.5) 



with 



tf'" " ^T^L-.*™ + ^b )r ' (M ' u) ' (A ' 6) 



For the / = amplitude T°(s,t,u) = 3A(s,t,u) +A(t,s,u) + A(u,t,s), with A(s,t,u) = 
T n + n -^. n o n o, the J = 0, 2 partial waves in large- A^c ChPT are found to be 

^ {S) = S wfff + 24^ t( 55L 2 + VL 3 )s 2 - 8(40L 2 + HL 3 )sml 



+10(58L 2 + 14L 3 - 3L 5 + QL 8 )m 4 n ] , (A.7) 
1 

120^/4 



= TTTT^TI ( 5L 2 + L 3 )(« " 4m^) 2 . (A.8) 



A.2 Unitarization of final state interactions 

Similar to meson- meson scattering amplitudes, it is possible to use the optical theorem to 
write down a series of unitarity relations for production amplitudes. Nonetheless, contrary 
to scattering, where in the physical region there are only absorptive cuts in the s-channel, 
in production processes one may also have unitarity cuts on the t- and u-channel within the 
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physical phase-space. Thus, strictly speaking, these channels are present in the unitarity 
relation. However, for the reasons explained in the main text, only s-channel unitarity is 
considered here. The unitarity relation taking into account tttt final state interactions in 
the s-channel is written as 

ImMj(s) = p(s)7j(s)*Mj(s) , (A.9) 

with p(s) = A(s, m 2 , mfy 1 / 2 / s = y/l — 4m 2 /s. In Eq. ([A.9Q , it is made explicit that 
the two-pion system in the final state of rf —> rjTnr is always 1 = (assuming isospin 
conservation). There are several ways to satisfy the former unitary relation. The simplest 
one is based on the if -matrix formalism where 

KA ( a \\ - Mj{s)\ tvee , . 

-Mj(s)| K -matri x - ~ . , ^7 71 , (A.10) 

1 - ip(s)Tj{s) | tr ee 

with M.j(s)\tree and Tj (s)|tree calculated at tree level and thus real. Through this formal- 
ism one incorporates the imaginary part of the pion loops. However, it misses the real part 
of the logarithms that appear, for instance, in the ChPT calculation at the loop level. One 
may think of completing the imaginary part ip(s) appearing in the If -matrix formalism 
with the full logarithm that shows up in the two-propagator Feynman integral 

16tt 2 B (s) = C- p{s) log ^f\±^- , (A.ll) 

p(s) - 1 

with ImBo(s) = p(s)/lQir for s > 4m 2 . This can be then incorporated in a solution of the 
unitarity relation of a form similar to the N/D method |29|, |3(| ^lj] : 

\a I M M-j{s) tree , » 10 \ 

Mj(s) N/D = , A.12 

l - 167rjB (s)7j J (s)|tree 

The two-propagator Feynman integral Bq(s) is nothing else but the g(s) function used in 
the N/D method |2^, ^0|, 0. Actually, the integral is ultraviolet divergent and has a 



local indetermination C (a SL (so) in the N/D analysis In a quantum field theory 

framework, this divergence is made finite by means of the corresponding coupling in the 
mesonic Lagrangian. 

B. 77-77' mixing 
B.l State mixing 

Let us consider the two-dimensional space of isoscalar pseudoscalar mesons. We collect 
the SU(3) octet and singlet fields in the doublet t?^ = (7/8,7/1). The quadratic term in the 
Lagrangian takes the form 

£ = \d^ T B Kd^ B - \i&MSb , (B.l) 
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with 



The mass matrix elements are defined as 

Mf = M| +AM| , 

Ml = M 2 + Mf +AMf , ( B -3) 
Ml = Ml +AMI , 

o 

where Mq denotes the U(1)a anomaly contribution to the r/i mass, Mf (i = 8, 1) the O(5 ) 
quark-mass contributions to the octet and singlet isoscalar masses, 

Ml= l -{AM 2 K -Ml) , M 2 =±(2M 2 K + M 2 ) , (B.4) 

o o o o o 

and Ml= -^p(M|- - M 2 ), with M\ and Mf the kaon and pion masses at O(5 ) in 
the combined chiral and l/Nc expansion. AMf , AM 2 and AM§j modify the lowest order 
values of the mass-matrix elements. 

To first order in 5s, 5\ and Ssi, the kinetic matrix K, can be diagonalised through the 
following field redefinition 8 : 



VB = z 1 / 2T -fj = z l / 2T - (J 8 J , zW-k-zW t = i 2 , 

z i /2= -«y 8 i/2 \ 

" I -Ssi/2 1 - ft/2 i 



(B.5) 



In the f\ basis the mass matrix takes the form 

M 2 = Z l/2 ■ M 2 ■ Z 1/2T , (B.6) 

where 

Mf = Mf (1 - 5s)- Ml 5 81 + AM 2 , 

M 2 = (M 2 + M 2 )(l - <5i)- Ml 5 8 i + AM 2 , (B.7) 
Mfi = Ml (1 - («S 8 + 50/2) - (M 2 + Mf + M 2 )<W2 + AM 8 2 X , 

to first order in AM 2 (and products 5 x AM 2 ). 

8 In general Z 1 ^ 2 ■ PC ■ Z 1 ^ 2 ^ — I2 but in the case of large- Nc ChPT the matrix Z 1 ^ 2 is real since chiral 
loops start at 0(S 2 ) and are not considered. 
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The physical mass eigenstates are obtained after diagonalising the matrix Ai 2 with an 
orthogonal transformation 



R T ■ M 2 D -R, fj = R T -r]p 

_ ( cos Op — sin 9p \ 
\ sin 6p cos Op I 



R 1 



(Bi 



One gets the relations 



and 



M 2 + M 2 , = Ml + M? 



Ml - Ml 

tanOp = — — 

Ml 



M 2 q , - M 2 



{M 2 -Ml) 2 +AMl , 



M 2 - M 2 



Ml 



Ml 



Ml 



(B. 



(B.10) 



M 2 - Mf M 2 - M 2 
Therefore, the matrix from the bare to the physical basis is given by r\B = {R-Z l l 2 ) T -ijp 



with 
R-Z 1 ' 2 



cosO P {l - 5 8 /2) + smO P 5 81 /2 -sin0 P (l - 6t/2) - cos0 P 5 81 /2\ 
smO P (l - S s /2) - cosO P 5 8 i/2 cos0 P (l - 5i/2) - smO P 5 81 /2 J 
In large- N c ChPT, 



8L 



P 



5 Ml 



and 



AM 2 



AM? 



Si = ^ M 2 +Ai , 
16L 8 



t>81 = — rr- Ms 



P 



81 ! 



P 

16Lg 
P 



Ml + Ml , , 



M\ 



Mt 



81 J 



AM, 



81 



2M 



81 



In RChT, 



(B.ll) 



(B.12) 



(B.13) 



and 



ICdCm M 2 
Ml P 



Si 



8c d c m M 2 
Ml f 2 ' 



^81 



8c d c m Ml 
Ml p 



AM 2 
AM 2 



8c 2 



8c 2 



M|/ 2 
MJ/2 



M 8 4 + Mg 4 ! 



M 4 + Ml 



MlM 2 K 



(B.14) 



(B.15) 



The low-energy constants of large- Nq ChPT are related to the RChT couplings (if higher- 
mass pseudoscalar resonance contributions are neglected) through L5 = CdC m /Mg and 
L 8 = c 2 m /2M 2 || . 
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B.2 Decay constants 

The decay constants in the rj-rj' system are defined as matrix elements of axial currents (at 
the quark level = <n M 7 5 -^<?) 

(0|A«(0)|P(p))=zV2/^ (a = 8,0;P = V,v') ■ (B-16) 
Each of the two mesons has both, octet and singlet components. Consequently, Eq. ( |B.1(]| ) 



defines four independent decay constants, fp. In the convention of Ref. [41| 



fn fn \ / /8 COS 8 -/ O sin0 O 

= I I I ) = { I • (B - 17) 

fv'fv'J V /s sin ^ 8 /o cos O 

In the chiral Lagrangian formalism, the physical masses and decay constants are ob- 
tained from the part of the effective action quadratic in the nonet fields with the correlator 
of two axial currents, the former from the location of the poles, the latter from their 
residues. Following the procedure outlined in Ref. the decay constants are given by 
/» = f[(F^)- l K] a p , where f = f„ = 92.2 MeV at 0{5 ) and F is the matrix which simul- 



taneously diagonalises the kinetic and mass matrices of the Lagrangian in Eq. (B.l) 



fC = fU 2 F , M 2 = F ] M 2 D F . (B.18) 



From Eqs. (|Blpl| ), we get F = R-jZ 1 / 2 )~ x and then f% = f\R\Z x l 2 )'%. The matrix 



F is exactly the same that relates the physical and the bare fields, rjp = R - (Z 1 / 2 ) 1 -np- 
To first order in 5s, <5i and 5 8 i, the fp are written as 



/»// = cos Ml + W - sm9 P 5 81 /2 , 
/°// = -sin0p(l + Si/2) + cos9 P 5 81 /2 
f*/f = sin0 P (l + 5 8 /2) + cos0 P <W2 , 
/£// = cos0 P (l + St/2) + sme P d 81 /2 , 



(B.19) 



'</' 

and the two basic decay constants and two angles as f 8 = f(l + 5 8 /2), /o = /(l + 5i/2), 
8 = Op + arctan((5gi/2), and 9$ = Op — arctan(<5gi/2), respectively. 
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